B decays into radially excited charmed mesons 
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It has been recently argued that some longstanding problems in semileptonic B decays can be 
solved provided the branching ratio for the B — > £)'(*' semileptonic decays are large enough. We 
have studied these decays in a constituent quark model which has been successful in describing 
semileptonic and non-leptonic B decays into orbitally excited charmed mesons. Our results do not 
confirm the hypothesis of large branching ratios for the B — > D'^*> semileptonic decays. In addition, 
we calculate the non-leptonic B — > D'n decays which can provide an independent test of the form 
factors involved in the B — ► D''*' reactions. 

PACS numbers: 12.39,Jh, 12.39.Pn, 13. 20. He, 14.40.Lb 



I. INTRODUCTION 

The determination of exclusive branching fractions of 
B -> X c l + vi decays is an essential part of the B-factory 
program to understand the dynamics of 6-quark semilep- 
tonic decays and therefore to determine the relevant 
Cabibbo-Kobayashi-Maskawa matrix elements, \V c b\ and 
\V u b\, which are essential parameters of the Standard 
Model. The methods used to determine the decay rates 
as a function of \V q b\ are very different for inclusive and 
exclusive decays, and the comparison of these comple- 
mentary determinations provides an important check on 
the results. 

Current measurements show a discrepancy between 
the inclusive rate B(B + -> X c l + vi) and the sum of the 
measured exclusive channels. At present, the disagree- 
ment between the inclusive rate B(B + —> X c l + v) = 
(10.92 ± 0.16)% and the sum of the measured exclusive 
rates is 



B(B + -> XJ+vi) - B(B + -> D^l+vi) 
- B{B+ -> D^irl+vi) = (1.45 ± 0.13)%. 



(1) 



The authors of Ref. [1] argue that this 1.45% discrep- 
ancy could be solved, or at least eased, by an unexpect- 
edly large B decay rate to the first radially excited D' and 
D'* states {B{B -> D'^l+ui) ~ 0(1%)). It is also stated 
in Rcf. [1] that a potentially large B{B + -> D'^l+vi) 
could help to solve the so called 1/2 vs 3/2 puzzle. This 
puzzle refers to the discrepancy between heavy quark 
symmetry based model predictions for B decays into IP 
D** states and observations. Those calculations predict 
that these decays should have a substantially smaller 
rate to the ji = -| doublet than to the ji° — f 
doublet (see for instance Refs. [2, 3]) whereas roughly 
equal branching ratios are found experimentally A large 
B(B + — > D'^l + vi) would result in an excess of the de- 
tected B — > D 1 / 2 ttIvi with respect to B — > D 3 / 2 ttIv due 
to the fact that the T(D' —> Di/ 2 tt) is much larger than 
the T(D' —> D 3 / 2 tt) because in the first case the pion is 



Resonance Mass (MeV) Width (MeV) 
L>(2550)° 2539.4 ±4.5 ±6.8 130 ± 12 ± 13 
D*(2600)° 2608.7 ±2.4 ±2.5 93 ± 6 ± 13 

TABLE I. Mass and total decay width of the D'^*' resonances 
as measured by the BaBar Collaboration [4]. 



emitted in S-wave while in the D 3 / 2 one is emitted in 
P-wave. 

Evidence for the D'W radial excitations, that would 
correspond to 2S* quark model states, have recently been 
found by the BaBar Collaboration [4]. Their results are 
summarized in Table I. The Z?(2550) meson has only been 
seen in the decay mode D*n and its helicity-anglc distri- 
bution turns out to be consistent with the prediction of 
a 2 1 S'o state. The £>*(2600) meson decays into Dir and 
D*tt final states, and its helicity-angle distribution is con- 
sistent with the meson being a J p = 1~ state. Moreover, 
its mass makes it the perfect candidate to be the spin 
partner of the D(2550) meson. 

Previous theoretical determinations of B{B + —} 
D'(*H + v) include an earlier calculation within the 
ISGW2 model, which incorporates constraints imposed 
by heavy quark symmetry, where a value of 0.06% was 
obtained [5]. This is in agreement with the heavy quark 
effective theory calculation of Suzuki et al. [6] which ob- 
tained B(B+ ^ D'^l+v) = 0.05%. A much larger value 
is obtained in Ref. [7] . There the authors use a relativis- 
tic quark model finding that the calculated branching ra- 
tio increases from 0.29%, for infinite heavy quark mass, 
to 0.40% when l/m q corrections are taken into account. 
Although the value is still far from the ~1% needed to 
explain the experimental discrepancy, the size of the cor- 
rections suggests that a complete calculation may further 
approach the 1% result. 

In this work we shall perform a full determination of 
B(B + -» D'^l + v) within the framework of the con- 
stituent quark model described in Rcf. [8]. The model 



has recently been applied to mesons containing heavy 
quarks obtaining a satisfactory description of many phys- 
ical observables: spectra [9, 10], strong decays and reac- 
tions [11, 12], and semileptonic and nonleptonic B and 
B a decays into orbitally excited charmed and charmed- 
strange mesons [13, 14]. We will use the factorization ap- 
proximation which gave a satisfactory explanation of the 
decays analyzed in Rcfs. [13, 14]. The form factors that 
parametrize the T(B —} D'lv) decay also appear, eval- 



L>(2550) as nj* 



2(T 



uated at q 



in the non-leptonic decay B — > D'it 



evaluated in factorization approximation. The latter de- 
cay, if experimentally accesible, could be used to extract 
information on the form factors near q 2 = 0. In this work 
we also evaluate the B° — > D i+ -k~ branching ratio and 
its branching fraction relative to the B° —t D + ir~ decay. 
The paper is organized as follows: In Sec. II, we de- 
scribe the constituent quark model predictions for the 
first radially excited 5-wave states paying special atten- 
tion to their strong decays. Sec. Ill is dedicated to ex- 
plain the theoretical framework through which we calcu- 
late the B{B+ -> D'^l+v) and B(B° -> D'+tt"). Our 
results are shown in Sec. IV. We summarize our conclu- 
sions in Sec. V. 



II. CONSTITUENT QUARK MODEL 
PREDICTIONS FOR 2S* STATES 



Channel T l 3p Bsp Q 

D*tt 131.90 99.87 

A> 0.18 0.13 

total 132.07 100 



TABLE II. Open-flavor strong decay widths, in MeV, and 
branchings, in %, of the D(2550) meson with quantum num- 
bers nJ p = 20". 



D*(2600) as nj p = 21" 


Channel 


F3 -Po 


&*, 


Dty 


10.84 


11.19 


D*it 


54.10 


55.83 


Dri 


11.86 


12.24 


D S K 


8.73 


9.01 


D*r/ 


9.65 


9.95 


DlTV 


0.28 


0.29 


D'nr 


1.44 


1.49 


DI-k 


0.01 


0.00 


total 


96.91 


100 



TABLE III. Open-flavor strong decay widths, in MeV, and 
branchings, in %, of the D*(2600) meson with quantum num- 
bers nj p —1\~. 



All the details on the constituent quark model we use 
have been described in Ref. [8] and we will only sketch 
here its main features. The model is based on the as- 
sumption that the constituent quark mass of the light 
quarks is due to the spontaneous chiral symmetry break- 
ing of the QCD Lagrangian. To restore the original sym- 
metry an interaction term, due to Goldstone-Boson ex- 
changes, appears between light quarks. This interaction 
is added to the perturbative One-Gluon Exchange (OGE) 
and the non-perturbative confining interactions. In the 
heavy quark sector, chiral symmetry is explicitly broken 
and Goldstone-boson exchanges do not appear. There- 
fore, the corresponding potential for the system stems 
from the nonrelativistic reduction of the OGE interaction 
and the confinement component. Explicit formulae and 
model parameters used herein can be find in Refs. [9, 13]. 

Concerning the new discovered mesons, D(2550) and 
£>*(2600), we predict masses of 2.70 GeV and 2.75 GeV 
which are larger than the experimental ones, 2.54 GeV 
and 2.61 GeV. While masses are large, the goodness of 
the meson wave functions has been tested through their 
strong decays. Those were evaluated in Ref. [11] in which 
a modification of the phenomcnological 3 Po decay model 
was proposed. The strength 7 of the decay interaction is 
scale-dependent being a function of the reduced mass of 
the quark-antiquark pair of the decaying meson. In this 
way a satisfactory global description of strong decays of 
mesons which belong to charmed, charmed-strange, hid- 
den charm and hidden bottom sectors was achieved [11]. 
Using physical masses for the 25 states, we obtained the 



results shown in Tables II and III. 

The D(2550) meson has been only seen in the decay 
mode D*7r, thus its possible spin-parity quantum num- 
bers up to J = 3 arc J p = 0~, 1 + , 2~ and 3 + . It is 
the lower in mass of the newly discovered mesons and 
within the possible assignments, the _ is the most plau- 
sible. Assuming it is a 25, J p = 0~ state, its total width, 
shown in Table II, predicted by the 3 Po model is in very 
good agreement with the experimental value given by the 
BaBar Collaboration. 

The D*(2600) meson decays into Dir and D*ir final 
states. Thus, its possible quantum numbers are J = 1~ , 
2 + and 3~. The hclicity-angle distribution of P*(2600) 
is found to be consistent with J p = 1~. Moreover, its 
mass makes it the perfect candidate to be the spin part- 
ner of the D(2550) meson. Again, the total decay width 
predicted by the 3 Po model, assuming it is a 25, J p = 1~ 
state, is in good agreement with the experimental data. 
Besides, we predict the ratio of branching fractions 



i3(P*(2600)° -)■ P H 
B{D* (2600)° -+D*' 



0.20 



(2) 



in reasonable agreement with experiment, 0.32 ± 0.02 ± 
0.09 [4]. 

In Tables II and III we also show the contribution of 
each channel to the total decay width of the D'(*> states. 
In both cases we obtain that the decay widths of these 
mesons into IP states are much smaller than the ones 
into 15 states. This is in contrast to Ref. [1] where the 
authors find plausible that the D'(*) decay rates to 15 



and IP charmed states may be comparable. That was 
used in Ref. [1] as a possible explanation of the so called 
1/2 versus 3/2 puzzle. However, we studied this issue 
within our model in Ref. [13] where we observed similar 
B semileptonic decay rates into the two j^ = 1/2+ and 
jq = 3/2 + doublets, being our results for the different 
channels in agreement with experiment. 



III. THE B + -> D'^l+vi AND B° ->• D' m ty DECAY 
WIDTH 

The presence of the heavy quark in the initial and fi- 
nal meson states in these decays considerably simplifies 

I 



their theoretical description. Let us start our analysis in 
the infinitely heavy quark limit, itlq — > oo. In this limit 
the heavy quark symmetry arises. This leads to a con- 
siderable reduction of the number of independent form 
factors which are necessary for the description of heavy- 
to-hcavy semileptonic decays. For example, in this limit 
only one form factor is necessary for the semileptonic B 
decay to S'-wave D mesons. It is important to note that 
the heavy quark symmetry requires that matrix elements 
between a B meson and an excited D meson should van- 
ish at zero recoil as a result of the orthogonality of the 
wave functions. 

As the D' and D'* states have 0~ and 1~ quantum 
numbers, the hadronic matrix elements for the semilep- 
tonic transition can be parameterized in terms of form 
factors as [15] 



(PV)|§ c (0) 7 ''(i-75) 1 M0)|fl(p)) 

(P / *(p')|fr c (0)7"(l - 75)*6(0)|g(p)) 
y/m B mD>* 



--h+(w)(v + v'Y + h-(w)(y - v'Y, 

--h v {w)e^etv' a vp 
- i[-h Al (w)(w + l)e*" + h A2 (w)(e* ■ v)v" + h A3 (w)(e* ■ v)v ,f *}, 



(3) 



where v(v') is the four velocity of the B(D'^*') meson, charmed meson and the form factors hi are dimensionlcss 
f oi23 _ _^ e n j g a polarization vector of the final vector functions of the product of four-velocities w = v ■ v'. The 

differential dT/dw decay widths are given by [15] 



^(B+ -+ D'l + H) = 9MM™k (w 2 1) 3/2 r 3 (1 + r) 2 G 2 H; 
dw 4<57T J 



^{B+ -> D'*l + n) = 9M}^JL (w 2 _ l)3/2 fu , + V 2 

dw 4871-3 



■ (w 2 - if' 2 (w + l)V 3 (l-r*) 2 



Aw 1 - 2wr* + r* 
w + 1 (1-r*) 2 



F 2 (w), 



(4) 



where Gf is the Fermi decay constant, \V c b\ the modulus a b — > c transition and r^*' = m D ,(,) /tub- As for G(w) 
of the Cabibbo-Kobayashi-Maskawa matrix element for and F(w), they are given in term of the form factors as 



1-r 
G(w) =h + (w) - j— h-(w), 



F 2 {w) ={ 2(1 -2wr* +r* 2 ) 



w — 1 



1 w + 1 



(1 -r*)h Al (w) + (w - l)(h Al (w) - h As (w) -r*h A . 2 {w)) 



x I (1 - r*f + -^-(1 - 2wr* + r* 2 ) 
W + 1 



(5) 



I 

Similarly to the semileptonic B + — > D^*'l + vi case we heavy quark symmetry predicts that all form factors are 
have that, in the limit of very large heavy quark masses, given in terms of just one Isgur-Wisc function £,{w). One 



has in that limit 



h+(w) = h v (w) = h Al (w) = h A3 (w) = £(to), (6) 



while 

h_(w) = h,A 2 {w) = 0. 
From these results one obtains in that limit 

G(w) =F(w) =£(to). 



(7) 



(8) 



Heavy quark symmetry at the point of zero recoil (w = 
1) implies £(1) = in the B — > £)'(*) case since the radial 
parts of the wave functions for the D'M and B mesons 
are orthogonal in the infinitely heavy quark mass limit. 
Thus, the value of £(w) near zero recoil comes entirely 
from corrections beyond that limit. This is different from 
the B -> £>(*) case where £(1) = 1 in that limit. 

The results of the above paragraph can also be ex- 
perimentally tested by means of the B° — ► D' + ir~ non- 
lcptonic decay because in the factorization approxima- 
tion the width of this process is related to the form fac- 
tors involved in the semilcptonic decay. 

Following the suggestion of Rcf. [16] we calculate the 
ratio B(B° -> D'+tt-)/B(B° -*■ D+7T - ) which in factor- 
ization approximation is given by 




where X(a, b 7 c) = (a + b — c) 2 
factor is related with h± via 



Aab, and the f+(q 2 ) form 



U(d 2 



TUB 

mp> 
m B 



m B 
raw 
m B 
m D i 



h+(u) 
h-(u). 



(10) 



Note q 2 and w are related through q 2 
2m B m DW u. If the B(B° -> D 



'+, 



- ln B f '""DO ~~ 

branching ratio 
is measured experimentally, it could be used to extract 

f?- 



information on f+^ D (0). 



IV. RESULTS 

In the calculation we use physical masses for the D'W 
mesons. In Fig. 1 we show the different form factors 
evaluated in our model for the B — > D'(*> transitions. For 
the sake of comparison we also show in a different panel 
the corresponding ones for the B — > D^ transitions. 
We see that, even for the actual heavy quark masses the 
relations, in Eqs. (6) and (7) are approximately satisfied 
for the B — >• £>(*' case over the whole w range. Deviations 
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FIG. 1. Upper panel: Form factors for the B + ->■ D^H+Vi 
transition evaluated in our model. Lower panel: The same 
for the B + — > D l( -*H + vi transition. 



are expected due to the finite heavy quark masses and 
the difference between mt and m c . For the B —> D'M 
decays those differences arc magnified by the fact that in 
the infinite heavy quark mass limit we have £(1) = in 
this case. 

The F(w) and G(w) factors are depicted in Fig. 2. 
Our results at maximum recoil, -F(u> max ) = 0.2 and 
G(u> max ) = 0.1, are compatible with the estimates in 
Ref. [1] 

F(«w) = 0.25 ± 0.15, G( Wmax ) = 0.15 ±0.1 (11) 

obtained adapting the light cone sum rule (LCSR) cal- 
culation of Ref. [17]. 

Integrating the differential decay width we obtain 
B{B+ -> D'l+vi) = 0.012% and B(B+ -> D'*l+v{) = 
0.097%. We thus have B(B+ -> D'Wj+i/j) = 0.109%. 
Although this branching ratio is larger than those found 
by Refs. [5, 6] it is still a factor ten smaller than the 
expectation in Rcf. [1]. By looking at the individual ra- 
tios we also find that B(B + — > D'l + vi) is smaller than 
B(B+ -> D'*l+n), in agreement with Refs. [5, 6]. 

Our results for G(w) and F(w) at zero recoil differ 
from the values obtained in Rcf. [7] using the rclativistic 
quark model. The discrepancy, clearly visible in Fig. 2, 
explains why our branching ratios are much smaller than 



0.3 



0.2- 



G(w) 

F(w) 

LCSR 

G(l)fromRef. [7] 

F(l)fromRef. [7] 



D u lv, 




FIG. 2. G(w) and F(w) factors for the B+ -> D l + v t and 
B + —¥ D'*l + vi transitions. We show with error bars the 
estimates for maximum recoil obtained in Ref. [1] adapting 
the light cone sum rule calculation of Ref. [17]. We also show 
the 1/m.Q corrected G(l) and F(l) values obtained in Ref. [7] 
using the relativistic quark model. 



the infinite heavy quark mass limit and neglects correc- 
tions on the inverse of the heavy quark masses. 

We have also evaluated the non-leptonic B(B° —> 
D' + ir~) branching fractions. The latter reaction may 
give additional information on the size of the form fac- 
tors involved in the semileptonic decay provided it can be 
measured in B-factorics or at LHCb in the near future. 
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the ones in Ref. [7]. The change is larger for F(l) which 
also explains why B{B+ ->• D'l+vi) > B(B+ ->• D'*l + n) 
in Ref. [7]. 

Concerning the non-lcptonic decay our theoretical re- 
sult is 



B(B° -> D'+n- 
B(B° 



D+n- 



= 0.011, 



(12) 



which combined with the experimental value for B(B° — > 
D+TT-) = (2.68 ± 0.13) x 10~ 3 gives for B(B° -> 
D' + ir~) ~ 3 x 10~ 5 , sufficiently large to be mesaured 
experimentally. 



V. CONCLUSIONS 

We have performed a theoretical calculation of the 
branching ratio for the B + — > D'^l + vi decays within 
the framework of the constituent quark model described 
in Ref. [8]. 

We find branching ratios much smaller than expecta- 
tions in Ref. [1]. Therefore, our results do not confirm 
that these contributions can explain the difference be- 
tween the inclusive B(B + — > X c l + v) rate and the various 
exclusive channels. 

Concerning the 1/2 versus 3/2 puzzle there is no need 
for a large branching ratio into the D'M states to solve it. 
In fact, it was already shown in Ref. [13] that our model 
predicts similar B semileptonic decay rates into the two 
jq = 1/2+ and jq = 3/2+ doublets, being our results for 
the different channels in agreement with experiment. To 
us, this apparent puzzle appears only when one works in 
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